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Abstract. Using single star models including the effects of shellular rotation 
with and without magnetic fields, we show that massive stars at solar metallicity 
with initial masses lower than about 20-25 Mq and with an initial rotation above 
~ 350 km s^^ likely reach the critical velocity during their Main-Scqucncc phase. 
This results from the efficient outwards transport of angular momentum by the 
meridional circulation. This could be a scenario for explaining the Be stars. 
After the Main-Sequence phase, single star in this mass range can again reach 
the critical limit when they are on a blue loop after a red supergiant phase 
(Hegcr & Langcr 1998). This might be a scenario for the formation of B[e] stars, 
however as discussed by Langer & Heger (1998), this scenario would predict a 
short B[e] phase (only some 10* years) with correspondingly small amounts of 
mass lost. 



1. Link between Be, B[e] supergiants stars and rotation 

A common feature of Be and B[e] supergiants is the non-sphericity of their 
circumstellar envelopes (see e.g. the review by Zickgraf 2000). More precisely, 
in both cases, disks are supposed to be present, likely disk of outflowing material. 
How do these disk form ? How long are their lifetimes ? Are they intermittent ? 
Are they Keplerian ? Many of these questions have been discussed in this 
conference and are still subject of lively debate. A point however which seems 
well accepted is the fact that the origin of an axisymmetri c wind structure such a s 
a disk might be connected to the fast rotation of the star ( Peluoessv et al. 2 000^ . 



If correct, this connection between fast rotation and the Be and B[e] phenomena 
leads to the question, when such fast rotation can be encountered in the course 
of the evolution of massive single stars ? In this paper, we give some elements of 
answer to that question based on models accounting for the effects of shellular 
rotation ((Zahn 1992 ). 



2. The critical velocity 

The critical angular velocity corresponds to the angular velocity at the equator 
of the star such that the centrifugal force exactly balances the gravity. The 
critical angular velocity r^crit,! in the frame of the Roche model for computing 
the gravity due to the deformed star, is given by 
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Figure 1. Le]t : Evolution of the total momentum of inertia of different 
stellar models at solar metallicity. The hook along the tracks corresponds to 
the end of the Main Sequence phase. Only the beginning of the track is shown 
for the 60 M© stellar model. .• The continuous line shows Evolution of 
the surface velocity of a 20 M© at solar metallicity with Uini = 300 km s^"'^ 
according to the models computed by Meynet & Maeder (2003). The long- 
dashed line corresponds to the critical velocity given by Eq. The short- 
dashed line shows the surface velocity the star would have, if solid body 
rotation is assumed. 



where i?pc is the polar radius when the surface rotates with the critical velocity. 

Looking at Eq. (P), one can be surprised that the stellar luminosity does 
not appear. Indeed, we could expect that, in addition to the centrifugal accel- 
eration, the radiative acceleration would help in balancing the gravity. When 
the stellar luminosity is sufficiently far from the Eddington limit (see below for 
a more precise statement), it has been shown by Glatzel (1998) and Maeder &: 
Meynet (2000) that radiative acceleration does not play any role. Physically, 
this comes from the fact that when the star reaches the critical limit at the 
equator, the effective gravity (gravity decreased by the effect of the centrifugal 
acceleration) becomes zero there and the radiative flux, responsible for the ra- 
diative acceleration, tends also toward zero due to the von Zeipel theorem (von 
Zeipel 1924; Maeder 1999). 

In contrast, when the stellar luminosity approaches the Eddington limit, 
the radiative acceleration becomes a dominant effect. Why such a difference 
between the case far from the Eddington limit and the case near the Eddington 
limit ? One could indeed argue that even near the Eddington limit, when the 
critical limit is approached, the radiative flux becomes zero at the equator. 
This is correct, but another mechanism comes into play here: the fact that the 
Eddington limit is modifled when the star is rotating. Let us first recall that 
the classical Eddington luminosity is given by the expression 

LEdd = 47rcGM/K, (2) 

where k is the total opacity, L the luminosity, M the mass of the star and the 
other symbols have their usual meanings. Now, when the star is rotating, two 
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Figure 2. Left : Evolution of the surface equatorial velocities at the surface 
of 60 models at solar metallicity with -Uini = 300 (short-dashed line) and 
500 km (long-dashed line). The continuous line shows the evolution of 
the equatorial critical velocity given by Eq. . Right : Same as the left part 
of the figure, for 60 at Z = 10~^ with Viui = 300 (short-dashed line) and 
800 km s^^ (long-dashed line). 



important differences appear: first the Eddington limit varies as a function of 
the colatitude 6, second, it is decreased when the rotational velocity increases. 
The Eddington limit modified by rotation is given by (Glatzel 1998; Maeder & 
Meynet 2000) 

LEdd(f^) = 47rcGM - /k(0), (3) 

where Ji = M/V{uj) is the average density of the star, V{lo) being the stellar 
volume when the star is rotating and to is the ratio of the angular velocity to 
the critical angular velocity given in Eq. (^. 

Thus, when the star is near the Eddington limit, rotation may sufficiently 
decrease L-^dd to make it equal to the actual luminosity of the star. In that case, 
we say the star has reached the OT-limit and strong mass loss ensues. 

Now at which velocity does this occur ? To obtain it, one has to find the 
value of J7 such that L = LEdd(^), where LEdd(^^) is given by Eq. ©. This is 
equivalent to find such that 

K T 02 
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where we have added the subscript "max" to indicate that the critical limit will 
be reached first at the position on the surface where the opacity is maximum. 
If the value of Q satisfying this equality is higher than f^crit,! (from now on 
called the classical limit), then only the classical limit is relevant since it will be 
reached first. To see if Eq. Q can be fulfilled with < Ocrit,!) let us use the 
definition of p and of lo to write 

'^-^V'(a;), (5) 
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with V'{ui) = 4 ^l^j ■ From the Roche model, it can be shown that the quantity 

^Lu'^V' (uj) increases from to 0.361 when uj varies from to 1 (see Maeder 
& Meynet 2000). It means that if Fmax is strictly inferior to 1 - 0.361 = 
0.639, Eq. © cannot be fulfilled. If Tmax = 0.639, Eq. © can be fulfilled 
with Q = r^crit,! and if Tmax is superior to 0.639, then values of < r^crit,! can 
satisfy Eq. Q). A new expression for the critical velocity, valid when the star is 
sufficiently near the Eddington limit, can be derived. It is given by 

f^crit,2 = (^) ^crit,l^/i^5^ (6) 

where Rd^) is the equatorial radius for a given value of the rotation parameter 
UJ. These expressions for the critical velocities are different from the expression 

f^crit = ^^crit,l(l — E) (7) 

used by some authors. Expression ((T)) is correct only if the surface is uniformly 
bright, which is not the case when the star is rotating fast. 

3. Evolution of the surface velocity: two simple extreme cases 

The evolution of the rotational velocity at the surface of stars depends mainly 
on three physical processes: 

• The efficiency of the angular momentum transport mechanisms in the in- 
terior, 

• The movement of expansion/contraction of the surface, 

• The mass loss. 

An extreme case of internal angular momentum transport is the one which 
imposes solid body rotation at each time in the course of the evolution of the 
star. A strong coupling is then realised between the contracting core and the 
expanding envelope. In that case, the angular velocity $7 is given by the ratio of 
the total angular momentum J and the total momentum of inertia of the star 
/. Eig. n shows the variation of / as a function of the growing stellar radius for 
a few stellar models at solar metallicity. In the case of the 9 Mq model, / varies 
as R" with a ~ 1. In the case of the 60 Mq model, a becomes negative. This 
results from the strong mass loss experienced by this star. 

Eor the 9 M©, since mass loss by stellar winds remains very modest during 
the Main Sequence phase, the total angular momentum is conserved during this 
phase. As a consequence, Q varies as the inverse of I and since I oc i?, oc 1 /R. 

Erom the previous section, we saw that when the star is sufficiently far 
from the Eddington limit, which is the case for the 9 Mq stellar model, then 
^crit,i oc R~^/'^ . Thus the critical angular velocity decreases more rapidly than 
the surface angular velocity when the star expands. Clearly this favors the 
reaching of the critical limit (see also Sackmann and Anand 1970; Langer 1998). 
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Figure 3. Left : Internal distribution of the angular velocity ri(r) as a 
function of the Lagrangian mass in solar units in a 15 Mq model, without 
magnetic fields, at various stages of the model evolution indicated by the 
central H-content Xc during the MS-phase. The initial velocity Uini — 300 
km s~^. Right : Same as the left figure but with magnetic fields. One notices 
the almost constant values of the angular velocity in the models. 

To illustrate this last point, let us consider the rotating track for the 20 Mq 
model (solar metallicity and initial rotational velocity of 300 km s~^) computed 
by Meynet & Maeder (2003) for obtaining values of the momentum of inertia 
during the evolution. From these values of / and also using the values for the 
actual total angular momentum, we deduce the surface velocity that the star 
would have in case of solid body rotation. We obtain the short dashed line in 
Fig. ^ Although the model is not self consistently computed (the evolutionary 
tracks were not computed imposing solid body rotation), it illustrates the fact 
that indeed when solid body rotation is achieved, the star may reach very easily 
the critical limit (here represented by the long-dashed line). 

Another extreme case is the case of no transport of angular momentum. 
Each stellar layer keeps its own angular momentum. The variation of is then 
simply governed by the local conservation of the angular momentum and, at 
the surface, oc 1/B?. When the radius increases, the surface angular velocity 
decreases more rapidly than the classical critical angular velocity, thus the star 
evolves away from the critical limit. 

Reality is likely in between the cases of solid body rotation and of local 
conservation of the angular momentum. Let us now see what are the predictions 
of more physical models. 

4. Evolution of the surface velocity in models with shellular rotation 

In the interior of stars, at least three mechanisms can transport angular momen- 
tum along a radial direction: 

• Convection:here we suppose that convective zones have solid body rotation 
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Figure 4. Evolution of the rotation velocities at the surface of 15 Mq models 

during the MS phase with and without magnetic fields. The initial velocity 
fini = 300 km s~^ in both cases. We see the much higher surface rotation 
when magnetic field is included. 



• Meridional circulation:this is the main mechanism for the transport of the 
angular momentum in radiative zones. 

• Shear turbulence:only the secular shear turbulence, occurring on thermal 

timescales, appears to be important. Dynamical shear, occurring on dy- 
namical timescales, appears in the advanced stages of the evolution of 
massive stars and only affects very locally the profile of J7 (Hirschi et al. 
2003). The efficiency of the vertical secular shear turbulence for trans- 
porting the angular momentum is in general much smaller than that of 
the meridional currents. 

Let us emphasize that the evolution of the angular velocity inside the star de- 
pends on the gradients of the chemical species and of the angular velocity, these 
gradients being themselves deduced from the variation of J7 as a function of the 
radius. Thus the problem has to be solved self-consistently, which is done in 
the computations presented here. More precisely, in the models discussed in 
this section, the effects of the centrifugal acceleration in the stellar structure 
equations are accounted for as explained in Kippenhahn & Thomas (1970) (see 
also Meynet and Maeder 1997). The equations describing the transport of the 
chemical species and angular momentum resulting from meridional circulation 
and shear turbulence are given in Zahn (1992) and Maeder & Zahn (1998) (sec 
details on the derivation of the angular momentum transport equation in Meynet 
8z Maeder 2005b). The expressions for the diffusion coefficients are taken from 
Talon & Zahn (1997), Maeder (1997). The effects of rotation on the mass loss 
rates is taken into account as explained in Maeder &; Meynet (2000). 
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Let us stress also that these models are able to account for many observa- 
tional constraints that non-rotating models cannot account for: they can repro- 
duce surface enrichments (Heger &: Langer 2000; Meynet &: Maeder 2000 ), the 
blue to red supergiant ratios at low metallicity ()Maeder Sz Mevnet 200lh . the 
variation with the metallicity of the Wolf-Ra vet popul ations and of the number 
ratios of type Ibc to type II supernovae (iMevnet &: Mae der 2005a). 

In the right part of Fig. ^ the evolution of the surface velocity for a 20 
Mq stellar model at solar metallicity with fini= 300 km is shown (see the 
continuous line) . Interestingly we note that the evolution of the surface velocity 
given by consistently taking into account the above transport mechanisms is not 
very far from the solid body rotation case, except at the very beginning and at 
the end of the Main Sequence phase. 

At the beginning, the differentially rotating model presents a decrease of 
the surface velocity, not shown by the solid body rotation model. This ini- 
tial decrease is due to the action of the meridional currents, which build up a 
gradient of inside the star, transporting angular momentum from the outer 
regions to the inner ones. This slows down the surface of the star. Then, in the 
interior shear turbulence becomes active and erodes the gradients built by the 
meridional circulation. Under the influence of these two counteracting effects 
the Q-profile converges toward an equilibrium configuration. This occurs on a 
very small timescale (a few percents of the Main-Sequence lifetime, see Denis- 
senkov et al. 1999). After this short phase, the variation of in the radiative 
zone continues to be shaped by shear turbulence, meridional circulation and the 
change of stellar structure (expansion/contraction of the stellar layers). 

At the end of the Main Sequence phase, when the star is older than about 
9 Myr, the surface velocity rapidly decreases. This is a consequence of the 
mass loss rate recipe we used in this computation (Vink et al. 2000; 2001) 
which shows important enhancement of the mass loss rates when some critical 
effective temperature are crossed (bistability limits, see the above references). In 
absence of such strong stellar winds, the surface velocity would increase during 
this phase. 

From this computation we can deduce the following results: first, during the 
Main-Sequence phase the transport mechanisms are efficient enough to maintain 
a relatively weak gradient of 17 inside the star (see the left part of Fig. 
The situation is thus not too far from the solid body rotation case. Let us 
however note that the gradients of although modest, are sufficient enough to 
drive chemical mixing. These models predict changes of the surface abundances 
during the Main-Sequence phase well in agreement with what is observed (Heger 
and Langer 2000; Meynet & Maeder 2000; Maeder and Meynet 2001). Second, 
this numerical example illustrates the importance of the mass loss in shaping 
the evolution of the surface velocity (see also the discussion below and Fig. I^J . 
Third, such models, starting with a higher initial velocity (typically above ~ 350 
km s~^) would easily reach the critical limit during the Main-Sequence phase 
(see the results shown is Meynet and Maeder 2005b). 

What does happen after the Main-Sequence phase ? The evolution speeds 
up and the variation of inside the star is mainly governed by the local con- 
servation of the angular momentum. We just saw in Sect. 3 above that, in 
situation where the radius grows up, this makes the surface velocity to evolve 
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Figure 5. Left : Evolution of rotating kinetic energy, of excess energy in 
the shear for a 15 M0 steUar model at solar metallicity, with Vi^i = 300 km s~^ 
without magnetic field. The energies are given as fraction of the gravitational 
energy. Right : The same as the left part but for model with rotation and 
magnetic field. The energy in the magnetic field is also indicated. 

away from the critical limit. Only when the star contracts, for instance when 
a blue loop occurs in the HR diagram, may the star reach the critical hmit 
( Heger &: Langer 19981 ) . This would be a possible scenario for the occurrence of 
B[e] stars, however as discussed by Langer & Heger (1998), this scenario would 
predict a short B[e] phase (only some 10^ years) with correspondingly small 
amounts of mass lost. 

Since such stars would have had their surface abundances changed by the 
deep outer convective zone appearing at the red supergiant phase, one expects 
that their surface abundance be highly enriched in CNO-processed material: as 
a numerical example, the N/C and N/0 ratios at the surface of a 9 Mq stellar 
model at solar metallicity, with an initial rotation of V[ni = 300 km s^^ are 
enhanced by factors equal to 2.5 and 2.0 respectively during the first crossing of 
the HR gap from the blue to the red. These two ratios become 8.5 and 5.3 on 
the blue loop after a red supergiant stage. The models without rotation would 
predict for these ratios the following values: first crossing, no enhancement for 
both ratios; on the blue loop, enhancement factors of 5 and 3.5 respectively 
((Mevnet Sz Maeder 2003^ . This numerical example shows that whatever the 
initial rotation velocity, models would predict some CNO-processed material at 
the surface if the star comes back from a red supergiant stage, secondly, rotation 
reinforces the surface enrichment on the blue loop with respect to non-rotating 
models. 

Let us now see what happens to more massive stars. In Fig. ^ the evolution 
of the surface velocity for various 60 Mq stellar models at two different metal- 
licities is shown. We see that at solar metallicity, the mass loss rates are so high 
that, even starting with an initial velocity of 500 km s~^, the star does not reach 
the critical limit. The situation is quite different at low Z, due to the metallic- 
ity dependence of the mass loss rates (here we used M{Z) = {Z/Zq,)"M{Zq) 
with a equal to 0.5 as devised by Kudritzki & Puis 2000). Indeed, at low Z 
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Figure 6. Time evolution of the surfaee helium content Ys in mass fraction, 
of the N/0, N/H and N/C in mass fraction for various models: The dotted 
line applies to model without rotation, the short-broken line to model with 
rotation (uini= 300 km s^^) but without magnetic fields, the continuous line 
to model with rotation (uini= 300 km s~^) and magnetic fields. 



little amounts of mass are removed by stellar winds, therefore little amounts of 
angular momentum. As a consequence, the angular momentum brought to the 
surface by the meridional circulation is not removed and accelerate the outer 
layers. Starting with vini = 300 km s~^, the star reaches the critical limit at 
the end of the Main-Sequence phase. Starting with an initial velocity of 800 km 
s~^, the reaching of the critical limit occurs at a much earlier time. 

For still higher initial masses, the stellar luminosity approaches the Edding- 
ton limit and, as explained above, the critical velocity becomes smaller than the 
value given by E g. (ITl). The star may rea ch the riF-limit and lose very high 
amounts of mass ()Maeder &: Mevnet 2000l ). Interestingly, in the HR diagram, 
the observed position of the de Jager or Humphreys-Davidson limit coincides 
with the position where this riF-limit would occur. This may be an indication 
that in the physics underlying this limit, both rotation and supra-Eddington 
luminosity play an important role. 



5. Evolution of the surface velocity in models with shellular rotation 
and magnetic field 

Spruit (2002) has proposed a dynamo mechanism operating in stellar radiative 
layers in differential rotation. This dynamo is based on the Tayler instabil- 
ity, which is the first one to occur in a radiative zone (Tayler 1973; Pitts & 
Tayler 1986). Even a very weak horizontal magnetic field is subject to Tayler 
instability, which then creates a vertical field component, which is wound up 
by differential rotation. As a result, the field lines become progressively closer 
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and denser and thus a strong horizontal field is created at the energy expense of 
differential rotation. 

In a first paper (|Maeder fc Mevnet 200.^ . we have shown that m a rotat- 
ing star a magnetic field can be created during MS evolution by the Spruit 
dynamo. We have examined the timescale for the field creation, its amplitude 
and the related diffusion coefficients. The clear result is that magnetic field an d 
its effects are quite important. In the second paper ( Maeder Mevnet 2004^ . 
a generalisation of the equations of the dynamo has been developed. The so- 
lutions fully ag ree with Spr uit's solution in the two limiting cases this author 
has considered ( Spruit 20021 ). i.e. "Case 0" when the /^-gradient dominates and 
"Case 1" when the T-gradient dominates with large non-adiabatic effects. Our 
more general solution encompasses all cases of fi- and T-gradients, as well as 
all cases from the fully adiabatic to non-adiabatic solutions. In a last paper 
(iMaeder fc Mevnet 200,^ . we examine the effects of the magnetic field created 
by Tayler-Spruit dynamo in differentially rotating stars. Magnetic fields of the 
order of a few 10^ G are present through most the stellar envelope, with the ex- 
ception of the outer layers. The diffusion coefficient for the transport of angular 
momentum is very large and it imposes nearly solid body rotation during the 
MS phase. This can be seen in Fig. |31 where are compared the evolutions of the 
angular velocity inside models with and without magnetic fields. 

The surface velocities resulting from these two models are shown in Fig. |31 
Except at the end of the Main Sequence phase, the model with magnetic field 
is strictly equivalent to the solid body rotation case. 

Does the model with magnetic field predict surface enrichments ? Shear 
turbulence in magnetic models is very weak due to the flatness of the internal 
profile. The excess of the energy in the shear in the magnetic model is only a 
few ten thousanths of the excess of the energy in the shear in the non-magnetic 
one (see Fig. On the other hand, solid body rotation drives meridional 
circulation currents which are much faster than usual and leads to much larger 
diffusion coefficients than the shear diffusivity and than the magnetic diffusivity 
for the chemical species. As a consequence, the surface enrichments obtained 
in the models with rotation and magnetic fields are higher than in models with 
rotation only (see Fig. ISJ. 



6. Conclusion 

Be stars might be the natural outcome of stars with initial rotational velocity in 
the upper tail of the initial velocity distribution. Depending on when the critical 
limit is reached, one expects more or less high surface enrichments. If the critical 
limit is reached very early during the Main-Sequence phase, no enrichment is 
expected, while if the critical limit is reached at the end of the Main-Sequence 
phase, high N/C and N/0 ratios are expected. At solar metallicity, for initial 
masses superior to about 50 Mq, mass loss rates prevent the stars to reach 
the critical limit. The lower initial value for reaching the critical limit is likely 
limited by variation of the distribution of the initial velocities. 

After the Main-Sequence phase, the variation of J7 inside the star is governed 
by the local conservation of the angular momentum. In phases during which 
the radius expands, this makes the surface velocity to evolve away from the 
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critical limit. In contracting phases, the reverse occurs. In the frame of single 
star evolution models, B[e] could be in a stage on a blue loop where the star 
contracted from a previous red supergiant phase. In that case, the surface is 
predicted to be enriched in CNO processed material. 

The effects of magnetic fields in this context remain to be studied. However, 
already at this stage, it appears that magnetic field will facilitate the reaching 
of the critical limit. 
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